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■ Abstract 

> 

^SJ ' In this study, we show that the energy eigenvalues and the eigenfunctions of the 

. Schrodinger equation for the power-law and the logarithmic potential can be easily 

obtained by using variation technique for special type wave functions. The results are 
in very good agreement with exact numerical results. 



1 Introduction 



(N 
O 

' A large number of important physical problems require solving the Schrodinger equation 

for spherical symmetric potential to determine the energy eigenvalues and the eigenfunc- 



tions. It is known that for very limited potentials, Schrodinger equation is exactly solvable. 
In general one has to resort to numerical techniques or approximation schemes, Most pop- 
ular approximation methods like expansion, WKB method, perturbation theory are 
widely used for this purpose. But some of these methods have-draw backs in application. 
Although some methods give simple relations for the eigenvalues, they give very complicate 
relations for the eigenfunction. The aim of present work is to give a simple way for finding 
both eigenvalues and the eigenfunctions of Schrodinger and Schrodinger-like equations for 
power-law and logarithmic potentials, which are very important in particle physics[l-10]. 
This paper is organized as follows. In the first section, the eigenvalues are obtained for 
the sgn{v)Ar" {v > —2) type potential by using variation technique for the special form 
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of eigenfunction. In section 2, the same technique is apphed to the logarithmic potential. 
In the last section, we give some concluding remarks. 



2 The solution of the Schrodinger equation for the power- 
law potentials 

The radial part of Schrodinger equation for sgn{v)Ar'^ {v > —2) type potential is written 
as 

gir) = Eg{r) (1) 



-— + sgn{v)Ar^ + ^^ 



Using substitution r = (;^)''+2p, Eq. (1) can be written in the following form 



-— + sgn{v)p^+ 



9{p) = ^9{p) 



(2) 



1 

where e = It is well known that Eq. (2) have exact analytical solutions for 

u = —1 and 1/ = 2 and the forms of the solutions are: 
for 1/ = —1 (Coulomb potential) 



gip)^p'+'Expi-xp)L''J+\2xp), 
ioT u = 2 (Harmonic potential) 

g{p) « p'+'Exp (- {xpf) Ln' {2{xpf). 



(3) 



(4) 



where L is Laguerre polynomials. Prom the similarity of these solutions, the following 
solution can be proposed for arbitrary v 



g{p) ^ p'^^Exp (- {xpf) Ln" {2{xp) 



(5) 



where x and d variation parameters and they can be obtained by minimizing of e in Eq. 
(2) with respect to these parameters. So; 

de 



dx 



= 0, 



(6) 
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Using Eqs. (2) and (5), for e, we get 

eni{x,d) = cx^ + (8) 



and from Eq. (6), we have x = """"^ ,thus, we have found e^i as the following 
form 



1/ 2 

e„Krf) = (t' + 2)(-) (-) , (9) 



where; 



^fc=o ^m=o «fc«r»^r(/c + m + 
ELoE;;^=oafc«m^(A; + m + ^) ' 



ELoEm=o«fe«mr(fc + m + ^) ' 



s = {2l + 1){21 + d+l) + {k + m-{k- mf)S, (12) 

where aj are the coefHcients of the generalized Laguerre polynomials. The parameter d can 
be obtained from Eq. (7), but unfortunately this equation can not be solved analytically. 
However, for n = 0, Z = case, the value of d which minimized of e can be found easily. 



In Fig. 1, we present the dependence d on v. Behavior of the curve is similar to + 2 
function. But, this function does not fit exactly to the curve. Therefore, a correction 
factor is necessary. According to the our assumption, the correction factor must be equal 
to 1 at u = —1 and v = 2. Thus, we chose the correction factor in the following form; 

w = {l + tpf (13) 
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where p 



— Ti ; — • iJ' 



a\v^+a2V+az ' 



io, we write 



d = ^/v + 2w 



(14) 



Fitting this equation to the curve in Fig. 1, t, ai, 02, as and h constants are obtained 
as, 0.2075, 0.1381, 1..05, 2.484 and 0.08104 respectively. So, when any values of v is given, 
d values corresponding to v can be easily calculated by using Eq. (14). In Tables 1 and 
2, we present results of our calculations for the eigenvalues at different u. For comparison 

in these tables, we also present existing numerical solutions in literature .In Table 2, the 

— V 

calculations done using eq. (9) are multiplied by 2 «+2 factor because of consistency with 
results of ref[ll]. 

When we study on linear potentials, Eq. (2) has the following form for S states 



where, z = p — e. The equation given above is known as Airy equation and has exact 
solution in terms of Airy functions given as the following form 



Where e„ are the zeros of the Airy function and as given in Table 3. For S states, our 
prediction for eigenfunctions has the following form from Eq. (5) 



^ + z g{z)=0 



(15) 



(16) 




(17) 
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In Figures 2 and 3, our wave function and the exact wave function arc presented to- 
gether for n = 0, Z = and n = 4, Z = states respectively. We see that the agreement 
between two solutions is excellent. In addition to the wave functions, in Table 3, eigen- 
values of the linear potential for S states (n = 0, 1,2,.. and I = 0) are given and compared 
with the exact results which are well-known zeros of the Airy function. Similarly, in Table 
4, eigenvalues of r°'^ are calculated and compared with exact numerical results. The cal- 
culated results are in good agreement with exact numerical results and better than given 
in ref[17]. 

3 The logarithmic potential case 

Let us consider V{r) = log(r) potential which is very important in particle physics. The 

radial part of Schrodinger equation for logarithmic potential is written as 

g{r) = Eg{r). (18) 

The function log(r) at = can be written in following form 

log(r)^iK-l]. (19) 

Substitute Eq. (15) into Eq. (14) and then apply the method presented in previous section 
to Eq. (14) for d we found the value 1.43203 at = 0.00001 (see Eq. (13)). Thus, the 
eigenvalues of Eq. (14) are written as 

E^l = ^-- (20) 

where is given as in Eq. (8) . The calculated results for eni are presented in Table 4. In 
this table we also give the results of numerical solutions. The eigenfunctions of logarithmic 
potential can be written from Eq. (5) as 

/ j\ 2i + l 

g{p)^p'+^Exp(-{xpY)Ln' {2{xpY) (21) 
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where, p = — It is obvious that the eigenvalues obtained for the logarithmic potential 
are in good agreement with the results obtained from numerical solution. Also, in order 
to show the validity of the wavefunctions given in Eq. (17) obtained for the logarithmic 
potential, first of all, some energy levels are chosen. Their wavefunctions are obtained 
by Numerov's Method. Finally, these results are given in Figures 4 and 5 together with 
predictions. 

4 Conclusion 

In this paper, we have calculated eigenvalues and eigenfunctions of power-law and loga- 
rithmic potentials by using variational techniques for special type wave functions. Our 
results are in good agreement with existing exact numerical ones. However, for the higher 
values of i/, there is some difference between two approaches. However this differences not 
larger than that one coming from 1/N expansion. Moreover, we obtained that present 
method predicts not only the eigenvalues as well as the eigenfunctions of given potentials. 
The correspondence of wavefunctions for lower energy levels is very good agreement with 
numerical results as expected the except case n and / arc large. The method used in this 
study present an easy way of calculating both eigenfunctions and eigenvalues of power-law 
and logarithmic potentials. In spite of its simple structure, the method is very practical 
and the results of the method are in good agreement with exact results. 
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Table. 1. Comparison of this study results for the ground state (n = 0, Z = 0) of various 

power-law i)ol cut ials. 



V{r) 


This Work 


Numerical[ll] 


V{r) 


This work 


Numerical[ll] 


_^-1.5 


-0.29703 


-0.29609 




3 


3 


_^-l.^S 


-0.22027 


-0.22029 




3.45110 


3.45056 




-0.25 


-0.25 


^4 


3.80241 


3.79967 




1 


1 


^5 


4.09626 


4.33801 


^0.15 


1.32798 


1.32795 


^6 


4.35243 


4.54690 


^U.S 


1.83352 


1.83339 


^7 


4.58158 


4.71772 


^0.75 


2.10829 


2.10814 


^8 


4.79013 


4.92220 


^1.5 


2.70816 


2.70809 


^10 


5.16092 
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Table. 2 Eigenvalues of — ^''^ and for different nl 



n 




This work 


numerical [1 1-15] 


n 


I 


This work 


numerical [11-15] 








-2.6859 


-2.686 








-1.2186 


-1.218 


1 





-2.2530 


-2.253 


1 





-0.4622 


-0.462 


2 





-2.0440 


-2.044 


2 





-0.2648 


-0.265 





1 


-2.3449 


-2.345 





1 


-0.5004 


-0.500 


1 


1 


-2.1006 


-2.101 


1 


1 


-0.2806 


-0.281 


2 


1 


-1.9504 


-1.951 


2 


1 


-0.1873 


-0.187 





2 


-2.1562 


-2.156 





2 


-0.2947 


-0.295 


1 


2 


-1.9900 


-1.990 


1 


2 


-0.1949 


-0.195 


2 


2 


-1.8749 


-1.875 


2 


2 


-0.1420 


-0.142 





3 


-2.0291 


-2.029 





3 


-0.2019 


-0.202 


1 


3 


-1.9049 


-1.905 


1 


3 


-0.1463 


-0.146 


2 


3 


-1.8124 




2 


3 


-0.1128 
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Table. 3 Eigenvalues of linear potential for different n and I = 0. 



n 


I 


This work 


numerical [16] 








2.33825 


2.33810 


1 





4.08918 


4.08795 


2 





5.52132 


5.52056 


3 





6.78614 


6.78671 


4 





7.94189 


7.94413 


5 





9.01859 


9.02265 
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Table.4 Eigenvalues of r°'^ for different nl together with exact values and other 



r(^soarch(-r rosuKs. wUli i")('rcou(agc- c-rrors 



n 


/ 


This work 


Numerical [17] 


Ref[17] 


% 








1.83352 


1.83339 


1.83375 


0.007 


1 





2.55152 


2.55065 


2.55142 


0.03 


2 





3.05177 


3.05118 


3.05224 


0.019 


3 





3.45197 


3.45213 


3.45341 


0.005 


4 





3.79233 


3.79336 


3.79482 


0.027 





1 


2.30056 


2.30050 


2.30073 


0.003 


1 


1 


2.85473 


2.85434 


2.85486 


0.014 


2 


1 


3.28666 


3.28583 


3.28659 


0.025 


3 


1 


3.64838 


3.64739 


3.64835 


0.027 


4 


1 


3.96361 


3.96268 


3.96382 


0.023 





2 


2.65760 


2.65756 


2.65775 


0.002 


1 


2 


3.12048 


3.12033 


3.12077 


0.005 


2 


2 


3.50296 


3.50245 


3.50309 


0.015 


3 


2 


3.83338 


3.83254 


3.83336 


0.022 


4 


2 


4.12686 


4.12581 


4.12678 


0.025 





3 


2.95448 


2.95445 


2.95461 


0.001 


1 


3 


3.35764 


3.35759 


3.35798 


0.001 


2 


3 


3.70299 


3.70270 


3.70327 


0.008 


3 


3 


4.00796 


4.00737 


4.00810 


0.015 


4 


3 


4.28282 


4.28196 


4.28283 


0.020 





4 


3.21236 


3.21233 


3.21247 


0.001 


1 


4 


3.57275 


3.57275 


3.57310 


0.000 


2 


4 


3.88913 


3.88898 


3.88950 


0.004 


3 


4 


4.17308 


4.17268 


4.17335 


0.010 


4 


4 


4.43196 


4.46131 


4.43164 


0.015 
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Table. 5 Eigenvalues of logarithmic Potential for different nl 



n 


I 


This work 


numerical [1 1-15] 


n 


I 


This work 


numerical [1 1-15] 








1.0445 


1.0443 


3 





2.5957 


2.5957 





1 


1.6412 


1.6430 


3 


1 


2.7465 


2.7440 





2 


2.0134 


2.0150 


3 


2 


2.8801 


2.8800 





3 


2.2842 


2.2860 


3 


3 


2.9996 


2.9990 


1 





1.8485 


1.8474 


3 


4 


3.1071 


3.1070 


1 


1 


2.1513 


2.1510 


4 





2.8293 


2.8299 


1 


2 


2.3875 


2.3880 


4 


1 


2.9498 


2.9480 


1 


3 


2.5798 


2.5810 


4 


2 


3.0592 


3.0600 


2 





2.2903 


2.2897 


4 


3 


3.1592 


3.1590 


2 


1 


2.4917 


2.1910 


4 


4 


3.2512 


3.2510 


2 


2 


2.6629 


2.6630 


6 





3.1770 


3.1791 


2 


3 


2.8106 




10 





3.6411 


3.6427 
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Figure captions 

Figure 1. The dependence of d on i/. 



Figure 2. Comparison of our wavefunction and Airy function for linear potential {n = 0, 
I = 0),Dot lines represent Airy function, 

Figure 3. Comparison of our wavefunction and Airy function for linear potential (n = 4, 
I = 0),Dot lines represent Airy function, 

Figure 4. Comparison of our wavefunction and corresponding numerical wavefunction 
for logarithmic potential 

(n = 0,1 = 0) 

Figure 5. Comparison of our wavefunction and corresponding numerical wavefunction 
for logarithmic potential 
(n = 4, Z = 4) 
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